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Motivated by a beautiful demonstration of the Faraday's and Lenz's law in which a small 
neodymium magnet falls slowly through a conducting non-ferromagnetic tube, we consider the 
dynamics of a magnet falling through a superconducting pipe. Unlike the case of normal conducting 
pipes, in which the magnet quickly reaches the terminal velocity, inside a superconducting tube the 
magnet falls freely. On the other hand, to enter the pipe the magnet must overcome a large elec- 
tromagnetic energy barrier. For sufficiently strong magnets, the barrier is so large that the magnet 
will not be able to penetrate it and will be suspended over the front edge. We calculate the work 
that must done to force the magnet to enter a superconducting tube. The calculations show that 
superconducting pipes are very efficient at screening magnetic fields. For example, the magnetic field 
of a dipole at the center of a short pipe of radius a and length L « a decays, in the axial direction, 
with a characteristic length £ ps 0.26a. The efficient screening of the magnetic field might be useful 
for shielding highly sensitive superconducting quantum interference devices, SQUIDs. Finally, the 
motion of the magnet through a superconducting pipe is compared and contrasted to the flow of 
ions through a trans-membrane channel. 

PACS numbers: 41.20.Gz, 74.25. Ha, 07.55.Nk 



I. INTRODUCTION 



There is a beautiful demonstration of Faraday's and 
Lenz's laws which became very popular as a result of 
an easy availability of powerful rare earth magnets Q, |t| . 
The demonstration consists of a long pipe made of a con- 
ducting, non-ferromagnetic material, such as copper or 
aluminum, and a neodymium magnet which is allowed 
to fall through it. One finds that the magnet takes a 
very long time to traverse the pipe. In fact for a tube 
of about 2m in length, the magnet takes almost 25s to 
finish the trip [9|! On the other hand, a non- magnetic 
object of the same dimensions falls through the pipe in 
less than Is. It is quite amazing to observe the falling 
magnet from the top aperture, the magnet appears to 
be moving through a very dense fluid. In reality, air pro- 
vides only a negligible resistance, and what actually slows 
the magnet is the force produced by the eddy currents 
induced in the pipe. This force is proportional to the 
velocity of the falling magnet. When the drag force be- 
comes equal to the magnet's weight, acceleration ceases 
and the fall continues at a constant terminal velocity. 
For strong rare earth magnets, the terminal velocity is 
reached very quickly. Perhaps surprisingly, in view of 
the complexity of the problem, it is actually possible to 
perform a fairly simple calculation which agrees quanti- 
tatively with the terminal velocity observed experimen- 
tally 0]. Curiously, the calculation also predicts that the 
terminal velocity should be proportional to the electrical 
resistivity of the pipe's material. This suggests that if the 
pipe is an ideal superconductor, the velocity of the falling 
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magnet should vanish. One can easily see, however, that 
this conclusion can not be right. Suppose that a magnetic 
dipole is created inside an infinite superconducting pipe. 
During the process of creation, the magnetic field inside 
the pipe is changing and the electric current is induced 
on its surface. The surface currents screen the magnet's 
field and prevent it from entering into the interior of the 
superconductor. In the case of an ideal superconduc- 
tor (inertialess electrons) considered in the bulk of this 
paper, the penetration length is zero and both the mag- 
netic and the electric fields are perfectly screened. By 
symmetry it is also clear that for an infinitely long pipe, 
the magnetic field produced by the induced currents is 
maximum precisely at the location of the magnet. Since 
the magnetic force on a dipole is proportional to the gra- 
dient of the field it must, therefore, vanish so that the 
magnet will fall without any resistance. The theory of 
reference is not applicable to the perfect conductors 
because it was explicitly constructed to treat normal met- 
als for which the magnetic permeability is very close to 
that of vacuum. Furthermore, the rate of decay of the 
induced currents in such metals is very fast, compared to 
the magnet fall velocity, allowing us to neglect the effects 
of self-induction 9] . Clearly, neither one of these condi- 
tions is met in the case of super or ideal conductors which 
dynamically screen magnetic field from their interior. As 
the resistivity of the pipe metal is decreased, there will 
be a crossover from the terminal velocity found in ref. [j| 
for normal pipes to the free fall velocity inside perfectly 
conducting pipes. 

Although a magnet "created" in the interior of an infi- 
nite superconducting pipe will fall freely under the action 
of the gravitational field, it takes work to bring (create) 
the magnet inside the pipe in the first place. This is so 
because the magnetic field lines, which for a free dipole 
spread throughout the space, must now be confined in 



a much more restricted volume. In this paper we will 
calculate the work that must be done to bring a mag- 
net into superconducting pipe of length L and radius a. 
Furthermore, the formalism developed here can be easily 
extended to study more general problems of screening of 
the magnetic field in cylindrical geometry 0, 3 , which 
are of particular interest for the development of reliable 
SQUID devices @. 



II. THE MODEL 

The model that we shall study is depicted in Fig. ^ 
A magnetic dipole, of moment m = mz, is brought from 
infinity and is inserted into cylindrical superconducting 
pipe of length L and radius a along the symmetry axis. 
This axis is taken to coincide with the z axis of the coor- 
dinate system. The radius of the pipe is assumed to be 
much larger than the wall thickness and will, therefore, 
be ignored. The Faraday's law of induction requires that 

where &(z, t) is the magnetic flux passing through a cross 
section of the pipe at position z and time t, and E is 
the local electric field. Since the tangential component 
of the electric field is continuous across the supercon- 
ductor/air interface, the right hand side of Eq. Q must 
vanish because no electric field can be present inside a 
perfect conductor. Thus, the flux passing through any 
cross section of the pipe must be constant in time. Fur- 
thermore, when the magnet is at infinity, the flux entering 
the pipe is zero, which then means that $(z,t) = at 
all future times as well. Vanishing flux is a direct conse- 
quence of the physics of cylindrical superconductors and 
must be used as a boundary condition for the solution 
of the Maxwell's equations. In Fig. 1 we have made an 
attempt to represent the magnetic field lines which are 
not allowed to go through the pipe because of the restric- 
tion on flux. Again we stress that this behavior is very 
different from what happens with normal metals |(| for 
which the electric field does not vanish and magnetic flux 
changes through different cross sections of the pipe. 

Inside the superconductor, magnetic field is zero and 
the continuity of the normal component of B requires 
that the normal component of the external field also van- 
ishes at the superconducting wall. The magnetic field 
lines must, therefore, be tangent to the pipe's surface. In 
general, however, vanishing of B„ at the interface is not 
sufficient to fully specify the boundary condition neces- 
sary for the existence of a unique solution to Maxwell 
equations. 
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FIG. 1: Schematics representation of a magnetic dipole de- 
scending into a superconducting tube. The flux lines are re- 
pelled from the tube. 

III. UNIQUENESS 

What are the boundary conditions which will make the 
field equation 

V x B = /i J (2) 

have a unique solution? Suppose that Eq. J5J allows for 
two distinct solutions Bi and B2 for the same value of 
the current density J. Since V • B = 0, we can always 
define a vector potential A such that B12 = V x Ai^- 
Now, consider the integral over all space 

J SB-SBd 3 r = J V x SA • V x SA d 3 r, (3) 

where <5B = B2 Bi and <5A = A2 — Ai The integrand 
on the right-hand-side can be cast into more convenient 
form using the identity (V x <5A) 2 = SA- V x V x SA + V- 
(SA x V x SA). Recalling that the current sources for the 
fields 1 and 2 are identical: = fio(^2 — Ji) = Mo S3 = 
V x <SB = VxVx SA, the integrand of Eq. © reduces to 
a perfect divergence. Using the divergence theorem, the 
volume integral can now be transformed into an integral 
over the bounding surfaces S, which in our case are the 
cylinder and the spherical shell of radius R = 00. 

J SB 2 d 3 r = J da n • (SA x V x SA). (4) 

Since the dipolar field decreases rapidly with distance, 
the contribution to the integral in Eq. Q coming from 
the spherical shell at R = 00 vanishes, and S reduces to 
the surface of the cylinder. Using the cyclical invariance 
of the integrand, we see that Eq. (J2J has a unique solution 
(<5B(r) = 0) either if n x SA = or n x <5B = 0, on the 
pipe surface. To have a well posed problem it is, there- 
fore, not sufficient to specify only the normal component 
of the magnetic field, instead the tangential components 
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of the field at the interface must be provided. In a cylin- 
drical geometry the boundary conditions posed in terms 
of the vector potential are particularly useful. Azimuthal 
symmetry around the z axis, restricts the vector potential 
to have only one non zero component in the direction, 
A = Atf, e#(r, z). The flux through a cross section of the 
pipe is then 

$(2;)— J B-nda = / Vx A-hda — <j> A- dr = 2na A ( f,(a, z) 

(5) 

Specification of the flux passing through the pipe is, 
therefore, equivalent to the specification of A$ and guar- 
antees that Eq. has a unique solution in the cylindrical 
geometry. 



IV. A PIPE OF INFINITE LENGTH 

For superconducting pipes of L = 00, Eq. J2J can be 
solved analytically, while for finite L only numerical so- 
lution is possible. We start, therefore, with the L = 00 
case. 

The dipole of moment m = mi is located inside the 
pipe on the axis of symmetry at z = r = 0. In the 
pipe's interior there are no free currents and the Am- 
pere's law, Eq. (J2J, reduces to V x B = outside the 
magnet. The magnetic field can then be written as a 
gradient of a scalar function B = — V</?. This defines 
the scalar magnetic potential <p, which also satisfies the 
Laplace equation since VB = 0. The magnetic potential 
ip = (fa + find is produced by the point dipole, 
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fio mz 



47r(r 2 



r 2\3/2 ' 



(6) 



and by the currents induced on the surface of the su- 
perconductor. Since Eq. JBJ is a solution of the Laplace 
equation, so must be find, 



199 n2 

- ^-r — ifind + O z ipind = 0, 

r or or 



(7) 



for r < a. We next note that find must be odd in z, 
free of singularities, and must vanish as z — > 00. Under 
these conditions, the solution of the Laplace equation Q 
can be written in terms of a Fourier integral involving 
modified Bessel function of the first kind, 



<pind{r,z) = — / A(k)I (kr)sm(kz)dk. (8) 
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According to the results of the previous section, the func- 
tion A{k) will be uniquely determined by the condition 
that &(z) — 0, which means that no magnetic field lines 
are lost to the wall. For infinite superconducting pipe 
this boundary condition, is, therefore, equivalent to the 
vanishing of the normal component of the magnetic field 
B n (a,z) = on the wall. This means that the normal 
component of the induced field must cancel exactly the 
field produced by the dipole, 
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47r(a 2 + z 2 ) 5 /2 



(9) 



at the air/superconductor interface. Using the standard 
identities, Eq. @ can now be written in terms of the 
modified Bessel function of the second kind K\(x), 
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d r (pind(r,z) \ r=a = 
k 2 K\{ka) sin(kz)dk 7 



(10) 



Combining expressions © and HH and using 
dla{x)/dx = I\{x), enables us to calculate the function 



A(k) = 



Ho m kK\(ka) 
47r h(ka) 



(11) 



The magnetostatic scalar potential inside an infinite pipe 
is then 



ip(r,z) 



Ho mz 



47r(r 2 + z 2 ) 3 / 2 



Ho m 
2tt 2 



+ °° dkksHkz)^4l (kr) 



(12) 
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The first term on the right hand side of Eq. (|12|l is 
the potential produced by the point dipole located at 
r = 0, z = 0, while the second term is the magnetostatic 
scalar potential produced by the electric currents induced 
on the superconducting surface. For large z, Eq. 112|) 
simplifies to 



(p(r, z) 



Ho m \ z \ 
2i:a 2 z 



Mhr/a) klM 
Jo(fci) 2 



(13) 



where ki is the first root of the Bessel function J\; 
Ji(fci) = with k\ k 3.831. The axial magnetic field 
of a dipole inside a superconducting pipe, 



a * r > z >~ 2na 2 Jo(fci) 2 



(14) 



is, therefore, strongly screened, with a characteristic 
length £ = a/k-y ~ 0.26a. 
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To confine a magnet inside a superconducting tube of 
small radius costs a lot of energy. The work necessary to 
achieve this can be calculated using a charging process, 



W 



f 1 1 
J dXm ■ B md (\m) = -m • B md (m), (15) 



in which the dipole is charged from to its final value m 
while the induced field Bj„d responds accordingly. Using 
Eq. O we find 



W = 



27r 2 a 3 



h(x) 



0.797- 



(16) 



The condition that the magnetic field lines must exit 
the pipe on the same side on which they entered (&{z) = 
0), implies that the magnet can not "probe" the full 
length of the pipe. Thus, if the pipe length is such that 
L > a, its field screening properties should be identical 
to those of a pipe of L = oo. With this observation in 
mind we are now ready to study superconducting pipes 
of finite length. 



V. PIPES OF FINITE LENGTHS 

Consider a superconducting pipe of finite length L 
placed along the z axis, whose center coincides with the 
origin at z = 0. Results of the previous sections sug- 
gests that if L > a, shortly after the dipole finds itself 
inside the tube, the magnetic field configuration should 
be identical to that inside an infinite pipe and the axial 
force should vanish. However, little can be said about 
what kind of forces act upon the dipole as it enters or 
exits the pipe. Finite length case should, therefore, be 
examined with some care. 

From sections §2 and §3 we recall that the magnetic 
flux at any cross section of a superconducting pipe must 
vanish. This information will be central to the practical 
aspects of the theory from now on. As the dipole ap- 
proaches the pipe from z — > +oo, surface currents are 
generated over the pipe's surface. It is convenient to 
imagine that the length of the pipe is subdivided into 
uniform rings, each carrying a circulating surface current 
density j(z). The vector potential produced by these cur- 
rents has only <fi component and its magnitude is given 
by the linear superposition, 



A ind (r,z) 



afi 



-L/2 



L/2 



K(r,z,z')j(z')dz', (17) 



where the kernel 

K(r,z,z')= I dfce- fc|z - z,| Ji(fcr)Ji(fca) (18) 



is obtained from the field produced by one thin ring 
0- The flux generated by the surface currents at z is 
&ind{z) = 2-KaA in d{a,z), see Eq. (JJjJ. For r = a the 



kernel can be evaluates explicitly in terms of the hyper- 
geometric function, 



K{a,z,z') 



a 2 2 .Fi[3/2,3/2;3;- 



T 



2\z-z'\ 



/|3 



(19) 



Condition that through each cross section of the pipe 
the net flux vanishes, $(z,f) = 0, leads to an integral 
equation for the surface currents, 



-L/2 



K{z,z')j{z')dz' 



-L/2 



2n [(Zr, 



13/2 : 



where z m is the coordinate of dipole. To determine j{z), 
Eq. i|2U|) is solved numerically by first discrctizing the 
integral and then performing a matrix inversion. Once 
the current distribution is known, the vector potential 
A in d is calculated using Eq. (fT7| . Knowing A in d, the 
induced magnetic field on the axis of symmetry r = 0, 
Bi„d = z/r d r (rAi n d)\r=o an d the magnetic force on the 
dipole F = zd z (m ■ B i „ d )| z=ZmiI , = o can be easily eval- 
uated. We note that forces arising from flux trapping 
in type-II superconductors are neglected in the present 
discussion. For this to be a good approximation, the 
superconductor must have either high critical state cur- 
rent density J c or be subjected to only small magnetic 
field [UOHg. We will discuss this more in conclusions. 



A. Checking the accuracy of the numerical 
procedure 

The accuracy of the numerical procedure used to solve 
the integral equation Eq. (|2C)|> can be judged by compar- 
ing it with the analytical solution for an infinite pipe, 
Eq. (|12|) . In Fig. [5] we compare along the axis of 
symmetry for an infinitely long pipe with the field cal- 
culated using a numerical integration of Eq. (|20|l for a 
pipe of L = 10a with a dipole located at z m = 0. The 
agreement is perfect and attests to the reliability of the 
numerical solution. Furthermore, since the axial mag- 
netic field is maximum at the position of the dipole, the 
force on it will vanish, in agreement with our previous 
discussion. 

In Fig. |21 we also show that asymptotically the mag- 
netic field is well approximated by Eq. (|14l) . Thus, in 
agreement with the previous discussion, the axial mag- 
netic field is screened exponentially even for pipes of finite 
length, as long as L > a. 



B. Edge effects for superconducting tubes of finite 
lengths 

Satisfied with the accuracy of the numerical procedure, 
we can now use it to study the edge effects associated 
with the finite length superconducting pipes. In Fig. |3| 
we plot the magnetic force felt by a magnet as it moves 
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FIG. 2: Comparison of the numerical (circles) and the ana- 
lytical (full line) calculations for the induced magnetic field. 
The analytical result is for L — oo, while the numerical solu- 
tion is for L = 10a. Nevertheless, there is a perfect agreement 
between the two. We also compare the analytical asymptotic 
form (dotted line) Eq. 114H of the axial magnetic field, with 
the result of numerical integration (squares) for L = 10a, 
showing a clear exponential decay of the axial field, even in- 
side a finite superconducting pipe. 



from infinity into the interior of a superconducting pipe 
of length L = 10a. Panel (a) shows a strong repulsive 
force near the pipe entrance which vanishes rapidly as 
the magnet penetrates into the pipe. In the case of small 
neodymium magnets (weight 6g) used in our demonstra- 
tions of the Faraday's and Lenz's laws [!j and a supercon- 
ducting pipe of radius a = 7.85 mm we find the repulsive 
force to be sufficient to support a weight of 1 kgl In panel 
(b) we plot the work necessary to bring a magnet from 
infinity to a point z. Clearly, there is a large electromag- 
netic energy barrier that the magnet must overcome to 
enter the pipe. Note that soon after the magnet crosses 
the front edge into the interior of the tube, the poten- 
tial stabilizes at a plateau the value of which is in per- 
fect agreement with the work found to be necessary to 
confine a dipole inside an infinite superconducting tube, 
Eq. I|lu|) . This agreement is, again, a consequence of a 
very efficient exponential screening of the magnetic field 
by the surface currents, even inside pipes of finite length. 

As the magnet approaches a superconducting tube, it 
experiences a magnetic force which opposes its motion. 
The magnitude of this force can be estimated using a 
simple scaling argument. Magnetic field produced by a 
dipole decays with distance as B magnet ~ m/s 3 . If the 
superconducting pipe is sufficiently narrow (L 3> a) and 
the magnet is not too close to the front edge, the mutual- 
inductance effects between the sections of the pipe can 
be neglected and the induced surface current density at 
distance s from the magnet will be j(s) ~ B magnet ~ 
m/s 3 . On the other hand, each one of these current 
loops will produce a magnetic field at the position of the 
dipole, dB pipe ^ magnet ~ j(s)ds/s 3 ~ mds/s 6 , and will 
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FIG. 3: Force acting on a magnetic dipole and the associated 
potential energy as it moves into and out of a superconducting 
pipe. Negative force is repulsive. Although, the force peaks 
look symmetric, in reality the force decays algebraically out- 
side the tube entrance, and exponentially into the tube. 



result in a net force 



am 



Sd+L 



ds/s 7 



(21) 



where Sd = z m — L/2 is the distance of the magnet to the 
front edge of the pipe and a is the pipe's polarizability 
factor. The scaling function is found to be f(x) = 1 — 
cc 6 /(l -I- a;) 6 . Asymptotically, the force F z behaves as 




a < \s d \ < L, 



(22) 



Thus, if I > a, the force on a magnet decays alge- 
braically outside the pipe and vanishes exponentially af- 
ter the magnet enters the pipe. 

The polarizability a can be calculated by considering 
the far asymptotic limit of Eq. I|22|). Under this con- 
dition, the magnetic field varies only slightly over the 
length of the pipe, and the induced current can be cal- 
culated analytically by solving Eq. I|2(J|I . We find, 



mL 

2lTZ 3 



(23) 



This current will produce a magnetic field at the position 
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FIG. 4: Asymptotic force felt by a dipole when it is far from a 
superconducting pipe. Points are the result of the numerical 
calculation and the solid lines are the estimates obtained using 
Eq. 12 1 1 . Panel (a) is for L = 10a and panel (b) is for L = 20a. 



of the dipole resulting in a repulsive force 

3{j,orn 2 a 2 L 



F 



(24) 



Comparing Eqs. I|24() and (12211 we find the polarizability 
of the superconducting pipe to be a — /ioa 2 /Sir. 

In Fig. 01a) we compare the asymptotic force calcu- 
lated using the numerical solution of the integral equation 
(|2U|l with the estimate obtained using Eq. (|?T|l for a tube 
of L = 10a. The agreement is very good for z m 3> L, but 
poorer for smaller values of the coordinate. The problem 
is that L = 10a is not sufficiently large to well satisfy 
the inequality o < z m <C £, as demanded by the first 
asymptotic region of l|22l) . For larger values of L, such as 
L = 20 of Fig. OJb), the agreement already is much bet- 
ter. The abrupt decay of Fs d near the entrance of the 
tube is a consequence of the finite cross section of the 
pipe, which becomes important for small values of Sd- 



VI. CONCLUSIONS 

We have examined the forces and the fields produced 
by a small magnet as it enters into a superconducting 
pipe. For an infinitely long pipe, magnetic field produced 



by a magnet inside a tube is exponentially screened by 
the surface currents. 

In the case of finite superconducting pipes we find that 
the magnetic field is also exponentially screened, as long 
as the length of the pipe is larger than its diameter. The 
exponential screening of the magnetic field is a conse- 
quence of the Faraday's law which leads to vanishing of 
the magnetic flux inside a superconducting pipe. The flux 
lines can not go through the pipe and at any cross section 
the number of lines going down is the same as the number 
of lines coming up. The electromagnetic cost of confining 
a dipole inside a superconducting tube is, therefore, very 
large, since the field lines must be strongly compressed 
to fit inside the pipe. We calculate that in order to insert 
a small neodymium magnet of 6g into a superconducting 
pipe of radius 8mm requires 10N of force. Once the mag- 
net is inside the pipe, however, its motion will continue 
unhindered. At the exit, the magnet will be ejected with 
the same force that was required to insert it into the pipe 
in the first place. 

It is curious to compare the magnet's motion through 
a superconducting pipe, with the flow of ions through an 
ion channel Ion channels are water filled holes, re- 
sponsible for a small potential gradient that exists across 
all biological membranes. Since the dielectric constant 
inside a pore is much larger than the dielectric constant 
of a phospholipid membrane, the normal component of 
the electric field at the channel/membrane interface is 
very small and the electric field lines are mostly confined 
to stay in the pore's interior. In this respect, the pore 
is very similar to a superconducting pipe, for which the 
normal component of the magnetic field also vanishes on 
the boundary. Nevertheless, one finds that contrary to 
what happens to a magnet inside a superconductor, the 
repulsive force on an ion does not vanish even when it 
is far inside the channel. Ions always feel a force which 
tries to expel them from the channel (except exactly at 
midpoint where, by symmetry, the force vanishes). 

The difference between superconducting pipes and ion 
channels is precisely due to the additional constraint of 
vanishing flux imposed by the Faraday's law on a per- 
fect conductor. The requirement that the flux through 
any cross section of a superconducting pipe must vanish, 
forces the system to be in a metastable state. If this con- 
dition is relaxed, energy of the magnet/superconductor 
system can be lowered. For example, suppose that the 
magnet is placed inside a normal conductor and only 
later the temperature is lowered until the superconduct- 
ing state is achieved. In this case the magnetic field con- 
figuration inside the pipe will be different from the one 
when the magnet is placed into an already superconduct- 
ing pipe. There will no longer be a restriction that the 
flux must vanish and, in fact, the magnetic field lines at 
the moment that the pipe turns superconducting will be 
partially expelled, leaving behind a finite flux. The en- 
ergy cost of confining a magnet under these conditions 
will, therefore, be significantly lower. Inside an ion chan- 
nel, behavior of the electric field is very similar to the 
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case of a superconductor without a vanishing flux con- 
straint. It is precisely this constraint that is responsible 
for the exponential screening of the magnetic field found 
for superconducting pipes of even finite length. If the 
vanishing flux constraint is relaxed (by say the process 
of turning the pipe superconducting after the magnet is 
placed inside), the magnetic field will no longer be ex- 
ponentially screened inside a finite pipe, and the magnet 
would encounter a repulsive barrier, similar to the one 
faced by an ion inside a trans-membrane channel, and a 
non-vanishing magnetic force [Tl| 

In this work we have neglected the hysterectic ( flux 
trapping) forces typically present in type-II superconduc- 
tors. For small cylindrical neodymium magnets weighing 
6g and radius r — 6.35 mm moving through a supercon- 
ducting pipe of a ~ 8 mm, the magnetic field at the pipe 



surface is on the order of B ~ 0.1 T. Considering a strong 
pinning condition with a high critical state current den- 
sity J c ~ 10 8 A/m 2 , we estimate 0,0,^3 the drag force 
on the magnet to be only a small fraction ~ 0.0001 of the 
force required to enter the superconducting pipe. Thus, 
our results should not be much affected even if more re- 
alistic type II superconductors are used. 
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